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ABSTRACT 


A fuel -minimized trajectory is generated, satisfying initial and final 
boundary conditions for the problem of a constantly thrusting space vehicle 
which is transferred from an initial elliptical orbit around Mars to a circular 
orbit of lower energy. The constant-thrust solution is amenable to electric 
ion-propelled spacecraft having ah initial elliptical or circular planetocentric 
orbit which must then be powered to a circular orbit of desired altitude. 
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Symbol 

<*> 


DEFINITION OF SYMBOLS 

Definition 


r 

u 

v 

T 

fi 

a 



I 

sp 


m 0 


m 


range angle (deg) 

radial distance of spacecraft to center of mass of planet (km) 
radial velocity (km/sec) 

velocity perpendicular to radial velocity (km/sec) 
thrust force (newton) 
steering control angle ( deg) 
radius of Mars (km) 

gravitational acceleration at surface of Earth 
( 9. 80665 x 10“ 3 km/ sec 2 ) 

gravitational acceleration at surface of Mars 
(3, 837 x 10“ 3 km/ sec 2 ) 


Mars gravitational constant 

( km 3 /sec 2 ) = g a 2 
m 

engine specific impulse 
(sec) 

initial mass (kg> 

constant time rate of 
change of mass (kg/sec) 


T 




t 


time (sec) 


Mars 



DEFINITION OF SYMBOLS (Concluded) 


Symbol 

m 

*2 

X3 

A 

h 


Definition 

mass at time t (kg) = m 0 - m(t - t 0 ) 

v 

first Lagrangian variable (sec) 
second Lagrangian variable (sec/km) 

third Lagrangian variable (sec 2 /km) 

fourth Lagrangian variable (sec 2 / km) 

semimajor axis of starting ellipse (km) 

coefficient of the second harmonic of Mar’s potential 
function J 2 = 0. 002001 



v 



USE OF LOW CONSTANT THRUST 
TO ACCOMPLISH ORBITAL TRANSFER 


SUMMARY 


A fuel-minimized trajectory is generated, satisfying initial and final 
boundary conditions for the problem of a constantly thrusting space vehicle 
which is transferred from an initial elliptical orbit around Mars to a circular 
orbit of lower energy. The constant-thrust solution is amenable to electric 
ion-propelled spacecraft having an initial elliptical or circular planetocentric 
orbit which must then be powered to a circular orbit of desired altitude. 


INTRODUCTION 

The Advanced Reconnaissance Electric Spacecraft (ARES) is a mission 
designed principally to map the neighboring planets of Venus, Mars, and 
Jupiter. To map the planets in a minimum length of time, high data transmission 
rates (on the order of 2 megabits/ see) are required, which, in turn, neces- 
sitate large quantities of electric power. The on-board power source is already 
required to accomplish the mission mapping objective once an orbit is estab- 
lished; hence, it appears practicable to use the available electric power source 
for main propulsion to achieve the desired orbit. For cartographic purposes, 
the most desirable orbits are circular, as opposed to elliptic, because the 
problem of changing altitude (and, hence, resolution) as a function of time is 
eliminated. In the ARES mission, once the planet is mapped from a high circular 
orbit, the altitude is decreased to a lower circular orbit so that a higher- reso*- 
lution map may be made for easy reference to the previous map. The optimal 
solution of the equations of motion for a minimum transfer time from one 
circular orbit to another is desired which includes low continuous thrust and 
the precessional effects of planetary oblateness. The results of such a solution 
for a typical problem are presented in this report. However, for purposes of 
generality the initial orbit is chosen as an ellipse, as if the spacecraft has been 
captured in an initial elliptical orbit and then powered to a circular orbit, 



DISCUSSION 


If 


The Differential Equation System 

The second order Newton differential equations of motion in plane polar 
coordinates are used, Included are the acceleration components from continuous 
thrust (T) along and perpendicular to the radius vector: 

r # * r<£ 2 + cos p - \x/v 7 ’ - | J 2 ^ (3 sin 2 ^ - i)J (ia) 


H sin p - (x/v 2 sin <£ cos <p j (lb) 

where the terms involving. J 2 within the brackets of ( i) account for the oblate- 
ness of Mars { 1) , 


I Introduction of the variables u » r and v - r<£ into ( 1) yields the following 

4 first order differential equation system: 
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‘V.j The variable (3 of set (2) is not associated with any particular differential 

| equation so it is removed by the definitions: 
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where X 8 , X 4 are functions of time. Recourse is made to the calculus of varia- 
tions to supply the Euler- Lagrange (E-L) differential equations for \ 8 and \ 4 
thus linking the steering angle p to the concept of minimal time. In this example, 
the solution satisfying the boundary conditions and the E-L equations yields a 
re lative sninimum* The E-L equations furnish the necessary criteria for a 
minimum to exist but are not sufficient to guarantee an absolute minimum* 
Another trajectory history might exist that uses less time than this example. 


The Appendix constitutes a brief exposition of the isoperimetric problem 
solved here. The integral to be extremalized specializes to the total time of 

M 

flight: I = J dt, f « i via equation (Al) . The constraining differential 
t 0 

equations are: 


Gi » 0 - X ss 0 , G« = U - — - — 
1 r r 9 3 r m 
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"7 j^ 1 - | Jj( j) (3 sin 2 <t> - 1)J 


as 0 


Go = r - u = 0. G. = v + — . — ^ — 

r m 


r 2 

+ 3 ^ sin (j> cos (p = 0 


The isoperimetric integrand function, equation (A5) , becomes 
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Application of the Euler-Lagrange equations (A-9) gives 

-9X 3 


sin <p cos (p 


9F sin <p cos <f> + 3X| (cob 2 <P - sin 2 
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The complete differential equation system is: 


0 =“ 
^ r 


r = u 
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\| = - 9 \ 2 sin <p cos i p + 3 X 4 *fy* a (cos 2 0 - sin 2 0 ) 
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Boundary Conditions 


The curves of Figure 1 illustrate to scale the elliptic starting orbit and 
the final circular orbit. Mars is located at a focus of the ellipse and the center 
of the circle. The radius of Mars is 3332. 15 km. The semimajor and semi- 
minor axes of the ellipse are 18 500 km and 16 000 km. The eccentricity is 
0. 502; the perihelion distance is 9212. 91 km; and the aphelion distance is 
27 787. 09 km. The radius of the circle is 13 332. 15 km, or at 10 000 km 
altitude above the surface of Mars. The initial spacecraft mass is 2950 kg, 

and the engine I is 3500 sec. The ionic thrust (T) is 0.4 newton. The 
sp 

constant time rate of change of mass (m.) is 1. 165 x 10’ 5 kg/sec. Initial and 
final conditions and derived quantities are given in Table 1. 


TABLE I. INITIAL AND FINAL CONDITIONS AND 
DERIVED QUANTITIES 


" ^Conditions 
Variables 

Initial Conditions 
and Derived Quantities 

Final Conditions 
and Derived Quantities 

<t> 

0 

78 orbits + 195° 

r 

9212.91 

13 332. 15 

u 

0 

0 

V 

2. 6355 

1. / 876 

t 

0 

64 days 4. 6 hours 

m 

2950. 0 

2885.37 

g 

501. 9 x 10" 6 

239. 7 x 10 -6 


* 
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Numerical Integration 


The procedure used was a powerful seventh order Runge-Kutta-type meth- 
od with automatic step size control developed by Dr. Erwin Fehlberg of the MSFC 
Computation Laboratory. A tolerance assuring eight-place accuracy throughout 
the history for all variables allowed a minimum step size of 227 sec and a 
maximum of 1355 sec. The first orbit required 115 time steps, or about 
78 000 sec (21.7 hr) ; the intervals are shown in Figure 1. The 78th orbit 
required 45 time steps each of 1085 sec, or about 47 000 sec (13. 1 hr) . The 
IBM 7094 computer was used with double precision ( 16-digit word length) , and 
machine time was 12 minutes for the trajectory history. 


Remarks 


Figures 1 depicts the trajectory history with motion in a counterclock- 
wise direction from the ellipse to the final circle. The first 50 orbits change 
the initial ellipse to a circular orbit about 14 700 km above the surface of Mars 
(Figs. 1A through IF) . The thrust vectors (small arrows) of Figures 1 A through 
IE never create more than a 90* angle with the longitudinal axis of a nonrotating 
spacecraft until the 51st orbit is reached. Orbits 51 through 78 shrink the radius 
of the circular orbit to the desired altitude of 10 000 km above the surface of 
Mars. The thrust vectors of Figures IF through 1H show continuous braking 
action, and they are nearly aligned with the tangent of the trajectory curve. 

The final boundary conditions are met about halfway through the 79th orbit. 

There is no evidence of precession as the ellipse is turned into the circle. 

Figures 2A through 2D give radius and steering angle as functions of 
time. Figure 2C shows the continuous braking effect of the thrust vector that 
oscillates about the angle (3 = 270° . 

The effects of oblateness were not significant for this example. The 
total time difference between two comparison runs for J 2 = 0 and J 2 = 0. 002001 
was slightly over 1000 sec, with about a one degree change in the range angle 

4>. 
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FIGURE 2A. RADIAL DISTANCE AND STEERING ANGLE p AS A FUNCTION OF TIME 
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FIGURE 2C. RADIAL DISTANCE AND STEERING ANGLE 0 AS A FUNCTION OF TIME 
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FIGURE 2D. RADIAL DISTANCE AND STEERING ANGLE 0 AS A FUNCTION OF TIME 


CONCLUSIONS 


The problem of changing the planetocentric orbit of a long-term, con- 
stantly low -thrusting spacecraft from an initial elliptic to a final circular orbit, 
so as to obtain a trajectory which requires a minimum of fuel is shown to be 
tractable. The effect of the oblateness of the planet Mars, which contributes 
perturbations to the equations of motion in the form of harmonics to the gravi- 
tational field representation, was shown to be negligible for the case considered. 
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APPENDIX 

The integral 

A 

1= J f(x 1 ,x 2 ,,..x K x 1 ,x 2 ,...x K , t)dt [2] 
to 


(A-l) 


is to be extremalized where the functions eligible for the extremalization must 
obey N ^ K differential equations as side conditions. The problem is defined 
as isoperimetric. 

A set of comparison functions X. = x^ + e r/ i (t) i = 1, 2, , . . K is selected 
with € an arbitrary parametric variable and 7 ^ (t) i = 1,2, . , . K arbitrary 
functions of t subject to the restrictions *?j(to) = ^(tj) =0, i = 1, 2 , , , . K, 

The integral becomes 


1(e) = / X K ,* t ,£ 2 ,... X K> t)dt, 

t Q 

and for I( e) to be an extremum for e = 0, it is necessary that I’(Q) - 0, Dif. 
ferentiation under the integral sign with respect to e produces 


ti 


I'(e) =/ 
*0 


3f 8X, , 9f 8X K , 8 f 8X, , , 3f 9X K 

ex, 8e ••• de -" ax.. 


‘K 


dt 


and for e = 0 this becomes 


I'(0) = / 


a L 


8 f 9f 

— r,, (t) + ... »I K 

K 


+ 


df 

dx< 


+ +~r rjj. (t) 

8x k K 


dt 


= 0 . 


(A -2) 
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N of the arbitrary functions r^( t) i = 1,2, . . . K are dependent because of the 
N differential equations of constraint 

G.(Xj,X 2, • • *X„, Xj,&2> • • • X__j t) = 0 j = 1,2, . . . N . 


Each of these N equations is differentiated with respect to e and multiplied by 
a Lagrangian variable to give 


V" 


de + '" + 9X r de 


+ !3!V| = 

9 x k 9e ^ ae 8 x k aej 


3 1, 2, i • » N 

and, for e = 0, this becomes 


‘0G. 0G. 0G. 0G. 

Mt) +..V+-T-J n K (t) + — r 1 ?L.(t)+... +— j Ot) = 

j L dx l dx K K ax, K 9x K 


) = 1,2, ... N 


(A-3) 


The equations of (A-3) do not disturb the numerical value of the integrand 
of (A -2) , so they may be added to it to form 


X'(0) = / 


N 9G.~1 T f N 

Jj *11 W + -** + oT" + Tj K 

j=l 3 dXl J L K j =1 J 


j (t) 0x K ?? K (t) 


9G." 

+ Jj M*) “T 3 ^l(t) + . . • 

L^i j=l 3 dx u 


77 k < t) > dt = 0 


(A -4) 


A new integrand function is defined to be 


F = f + ^ \.(t) G ' (A-5) 

j=l 1 * 


which when differentiated partially with respect to x. and x i - 1 , 2, . . . K 
gives 



The right side of these equations is identical to the bracketed terms of (A-4) 
so that 


I’(0) 


/•VeF 

" { V 8x ‘ 


TJi<t> + ... + 7J K (t) 

K 


. 8F . . . SF • 

+ — n l(t) +... + — rj (t)\ 
dx t 8 x k k ) 


dt = 0 


(A-6) 


Integration by parts on the right half of the integrand of (A-6) is performed 
using the formula 
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U i “ 


9F 

dx. 


dv A = rj^t) dt 


du. = v. = Tf. ( t) 

1 W./ 1 1 


i =1,3,... K 


u.v. 
i i 


OF m 

— T?.(t) 

8x. 1 


it 


j v du = / i, t» 


d / 0F 


dt 


0x. 


dt . 


Equation (A-6) becomes 
T(0) = 


“T“ ^l(t) + . . . + 

ax, ax 


k i 


T] K (t) 

aX K K J 


t 

- I 

to 


8F 

0Xj 


_d_ 

dt 



T)i (t) + . . . 


8F 

0X 


K 


_d_ 

dt 


8F 


0X 


K, 


r] K (t) > dt = 0 . 


(A -7) 


The first term of (A -7) is zero because of the restrictions 77 . (t 0 ) = i?.(tj) = 0 
i =1,2,... K so that (A-7) becomes 


I’(0) = 


I 

to 


9F 

0Xj 


_d 

dt 



77 1 ( t) + . . . -f 


8F 


0X 


K 


_d 

dt 



r? K ( t) >dt = a 


( A-8) 
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The Lagrangian variables A^(t) j =1,2, . . , N can be algebraically 

assigned to Insure that the coefficients of the first N terms of (A-8) will 
vanish. (The K - N terms can be unconditionally brought to zero because 
77^ t) i = N + 1, N + 2, , . . K are independent and arbitrary. ) The Euler- 

Lagrange differential equations of (A -8) 



(A 


provide the differential equations for the Lagrangian variables. 


